Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1140 1713 (N42W34 L=151)   271     69        2 (1107/1697-1180/1726) -31.0 W31-W36 N39-N43
      2 1570 1664 (N24W62 L=179)   785    189        0 (1489/1620-1635/1693) -41.6 W55-W66 N19-N27
      3 1455 1550 (N22W44 L=160)  1425    422        3 (1293/1508-1641/1610) -35.5 W34-W57 N14-N29
      4 1283 1331 (N15W23 L=140)   279    103        0 (1232/1321-1341/1343) -26.6 W20-W27 N13-N17
      5  975 1280 (N21W04 L=120)   282     85        0 ( 932/1264-1018/1294) -31.0 W01-W06 N19-N22
      6  425 1157 (N28E35 L= 81)   960    249        2 ( 321/1101- 523/1216) -29.4 E44-E27 N24-N33
      7  751 1071 (N15E14 L=102)   802    351        2 ( 621/1035- 900/1116) -19.4 E22-E05 N11-N18
      8 1754  869 (S24W43 L=160)   299     73        0 (1722/ 850-1791/ 890)  -8.1 W40-W47 S25-S23
      9 1453  687 (S26W17 L=133)  1362    359        5 (1379/ 527-1488/ 851)  56.3 W08-W21 S34-S16
     10  604  371 (S22E49 L= 68)   654    176        2 ( 533/ 349- 684/ 400) -31.0 E53-E42 S24-S18
     11  453  357 (S20E64 L= 52)  1043     78        1 ( 403/ 337- 498/ 379) -10.3 E71-E59 S22-S18
  

    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