Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  511 1718 (N23E58 L=196)  1900    288        7 ( 451/1632- 554/1808) -39.1 E69-E49 N16-N29
      2 1678 1482 (N39W34 L=287)   632     80        2 (1648/1448-1709/1521) -51.3 W30-W36 N37-N42
      3 1818 1341 (N35W49 L=303)  1121    176        0 (1735/1285-1884/1383)   0.0 W39-W58 N33-N36
      4  171 1324 (S07E71 L=183)  2073    134        4 ( 120/1265- 235/1376)   6.7 E79-E60 S10-S03
      5  256 1232 (S11E57 L=197)  1599    256        3 ( 224/1101- 299/1326) -81.3 E60-E52 S18-S05
      6  172 1231 (S12E67 L=187)  1116    137        3 ( 144/1176- 205/1294)  54.5 E70-E63 S15-S08
      7  494 1174 (S11E35 L=218)   278     73        0 ( 470/1136- 511/1209) -36.9 E37-E33 S12-S08
      8  797 1042 (S12E13 L=240)   178     93        0 ( 778/1011- 829/1086) -31.0 E15-E10 S12-S09
      9 1385  680 (S16W30 L=284)   212     86        0 (1360/ 636-1416/ 720)  90.0 W29-W31 S19-S13
     10 1526  612 (S16W42 L=296)   351     66        1 (1501/ 581-1545/ 652) -31.0 W40-W45 S17-S13
     11  856  432 (S46W07 L=261)  1986    436        4 ( 684/ 357-1042/ 509)   5.4 E08-W24 S49-S43
  

    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