Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  487 1674 (N62E31 L= 39)   601     41        1 ( 467/1649- 511/1697)  26.6 E36-E27 N59-N64
      2  375 1654 (N64E56 L= 15)   641     68        1 ( 347/1627- 408/1670)  -6.8 E70-E45 N61-N66
      3 1289 1534 (N26W31 L=102)   266     65        0 (1266/1504-1314/1563) -90.0 W30-W32 N23-N28
      4  864 1515 (N37W04 L= 74)   481    135        0 ( 806/1485- 929/1537) -45.0 W00-W07 N33-N39
      5  472 1479 (N46E29 L= 41)   917    161        0 ( 398/1448- 549/1500)  -3.6 E37-E21 N43-N47
      6  438 1187 (N27E34 L= 37)   623    168        0 ( 409/1143- 499/1253) -53.1 E36-E29 N23-N31
      7  185 1125 (N29E61 L=  9)   641    105        0 ( 163/1074- 204/1174) -49.4 E64-E58 N25-N32
      8  789 1131 (N16E11 L= 60)   194     73        0 ( 778/1083- 807/1164)  53.1 E12-E09 N13-N17
      9  159 1064 (N25E65 L=  6)   447     36        0 ( 145/1034- 172/1092) -63.4 E67-E62 N23-N27
     10  604 1037 (N15E24 L= 46)   616    149        3 ( 554/ 975- 650/1079) -81.9 E26-E21 N11-N18
     11 1523 1009 (S09W29 L= 99)   223     45        1 (1503/ 981-1544/1039)  56.3 W26-W30 S10-S07
     12 1501  741 (S24W22 L= 92)   838    208        2 (1465/ 672-1557/ 844) -90.0 W19-W25 S29-S17
     13 1695  686 (S34W37 L=108)   581     78        2 (1676/ 645-1721/ 734) -74.1 W35-W40 S37-S30
     14 1275  630 (S25W05 L= 75)   455    153        0 (1192/ 611-1341/ 646) -18.4 W00-W09 S26-S22
     15  622  312 (S29E51 L= 19)  2558    299        2 ( 558/ 197- 733/ 420)  17.2 E77-E35 S36-S23
     16 1449  340 (S53W12 L= 82)   637     98        0 (1398/ 324-1492/ 362) -21.8 W06-W16 S55-S50
     17 1315  233 (S58E06 L= 65)   659     48        1 (1288/ 211-1346/ 260) -60.3 E09-E02 S61-S53
     18 1221  178 (S60E19 L= 51)   530     53        0 (1201/ 158-1247/ 198)  20.6 E23-E14 S62-S56
  

    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