Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  949 1818 (N52E29 L=199)   433     45        1 ( 928/1795- 976/1841) -15.9 E32-E25 N49-N54
      2 1020 1723 (N45E16 L=212)  1158    133        3 ( 963/1689-1091/1764)   5.2 E21-E09 N42-N48
      3 1024 1607 (N36E11 L=216)   435     97        1 ( 985/1587-1075/1631)   9.5 E13-E07 N34-N37
      4  324 1589 (N22E71 L=157)   706     61        0 ( 300/1560- 344/1616)  -6.7 E82-E65 N20-N24
      5  274 1418 (N11E63 L=164)  1281    174        5 ( 240/1364- 297/1503) -66.0 E69-E58 N07-N16
      6  402 1471 (N16E52 L=175)   369     69        0 ( 377/1449- 436/1501) -14.0 E56-E48 N15-N17
      7 1305 1185 (N13W14 L=242)   858    190        3 (1235/1138-1366/1229)  45.0 W10-W18 N09-N17
      8  725 1189 (N04E21 L=207)   345    106        0 ( 687/1166- 773/1227)  -9.5 E23-E17 N03-N06
      9 1256 1113 (N08W12 L=240)   916    286        3 (1180/1044-1349/1220) -22.6 W06-W18 N04-N13
     10 1022 1129 (N06E02 L=226)   480    125        0 (1006/1066-1043/1199) -66.8 E03-W00 N02-N09
     11  372 1027 (S11E43 L=185)   762    246        0 ( 355/ 910- 404/1133)  90.0 E44-E41 S17-S03
     12 1542 1014 (N07W32 L=260)   262    101        1 (1501/ 975-1565/1038) -21.8 W29-W34 N05-N08
     13 1268  778 (S11W19 L=247)   971    348        0 (1196/ 651-1355/ 930)  81.0 W16-W21 S19-N00
     14  286  863 (S22E50 L=178)  1004    184        0 ( 207/ 827- 369/ 887)  14.7 E60-E41 S25-S19
     15 1643  615 (S15W50 L=278)   341     96        0 (1598/ 588-1673/ 637)  -7.1 W45-W53 S15-S13
     16 1472  424 (S30W46 L=274)   468    116        0 (1425/ 402-1529/ 454)  39.8 W42-W49 S31-S26
     17  605  361 (S53E24 L=203)  1862    364        2 ( 449/ 282- 718/ 449) -15.4 E41-E11 S58-S48
     18  954  234 (S56W15 L=243)   884    119        3 ( 907/ 196-1004/ 281) -18.4 W08-W23 S59-S52
  

    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