Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1577 1530 (N43W32 L=302)   344     67        0 (1552/1509-1603/1546)  33.7 W29-W35 N40-N44
      2 1329 1210 (N16W15 L=285)  3995    889        5 (1023/ 842-1589/1526) -32.0 E10-W39 S01-N29
      3  290 1343 (N05E59 L=211)   552     91        1 ( 250/1325- 322/1370) -45.0 E62-E55 N03-N06
      4  192 1047 (S15E61 L=209)   785    177        0 ( 164/ 999- 218/1077)  68.2 E64-E58 S17-S12
      5  126 1012 (S18E69 L=200)  1280    120        0 ( 108/ 964- 154/1054)  63.4 E72-E65 S20-S14
      6  235  985 (S18E55 L=214)   479     98        1 ( 216/ 939- 255/1031)  80.5 E56-E53 S20-S14
      7 1307  771 (S11W22 L=292)   735    180        1 (1206/ 730-1394/ 800)  32.5 W16-W27 S15-S07
      8  527  782 (S25E28 L=242)   348     81        0 ( 482/ 770- 569/ 796)   8.1 E31-E24 S25-S23
      9 1036  704 (S20W07 L=277)   737    166        0 ( 965/ 683-1132/ 742)   5.2 W02-W13 S21-S18
     10 1081  432 (S37W18 L=287)   251     36        0 (1064/ 411-1097/ 459) -76.0 W16-W19 S38-S34
     11 1316  312 (S40W42 L=311)  3008    504        1 (1087/ 283-1484/ 364)  12.3 W23-W56 S44-S33
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y