Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1415 1824 (N48W72 L=328)   766     63        0 (1395/1794-1443/1850)   2.3 W62-W87 N46-N49
      2 1528 1708 (N38W64 L=320)   674    135        0 (1487/1679-1574/1731) -51.3 W56-W68 N34-N39
      3 1535 1646 (N34W58 L=314)   512    109        0 (1498/1608-1575/1679) -80.5 W56-W59 N30-N36
      4  696 1517 (N43E16 L=240)   387     89        1 ( 658/1498- 752/1534) -20.6 E19-E11 N41-N44
      5 1558 1443 (N21W46 L=302)   312     79        0 (1533/1413-1590/1463)   8.1 W42-W49 N19-N22
      6 1751 1376 (N13W61 L=317)   570    117        1 (1725/1324-1780/1430)  82.9 W58-W62 N08-N16
      7 1261 1286 (N18W20 L=276)  1172    389        1 (1079/1229-1422/1331) -31.0 W08-W28 N11-N23
      8  467 1306 (N32E36 L=220)   352    114        0 ( 419/1280- 521/1321)  -6.3 E40-E31 N30-N33
      9  664 1226 (N24E20 L=236)   210     71        0 ( 636/1214- 709/1246) -21.8 E22-E17 N23-N25
     10  216 1202 (N28E62 L=194)   900     87        3 ( 186/1163- 255/1231)  33.7 E67-E57 N25-N29
     11  362 1128 (N22E45 L=211)   290     88        0 ( 329/1117- 418/1148) -12.5 E48-E39 N20-N23
     12  210 1101 (N21E62 L=195)  1173     83        3 ( 178/1063- 256/1138) -16.7 E66-E55 N18-N23
     13 1717  996 (S08W46 L=302)   660    126        1 (1692/ 937-1741/1052)  74.1 W42-W47 S11-S04
     14  149  898 (N09E72 L=184)  1608    129        4 ( 117/ 848- 179/ 950)  36.9 E80-E66 N05-N11
     15 1591  413 (S47W38 L=294)   718     97        2 (1553/ 384-1631/ 445)  23.2 W32-W42 S49-S44
     16  599  333 (S32E49 L=208)   640     89        2 ( 570/ 296- 622/ 377)  26.6 E55-E43 S35-S28
     17 1290  253 (S54W06 L=262)   452     55        0 (1268/ 234-1321/ 275) -26.6 W03-W09 S56-S51
     18  965  205 (S50E27 L=229)   408     21        1 ( 945/ 186- 988/ 222)   0.0 E29-E24 S52-S48
  

    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