Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  559 1560 (N28E35 L=225)   556    133        1 ( 539/1500- 577/1618) -76.0 E37-E33 N24-N32
      2 1690 1440 (N26W46 L=307)  1575    312        0 (1586/1300-1780/1533) -49.4 W38-W52 N17-N31
      3  443 1477 (N22E43 L=218)   447    126        0 ( 395/1451- 482/1497) -14.0 E46-E38 N20-N23
      4 1375 1419 (N22W21 L=282)   258     36        1 (1350/1398-1403/1437)  26.6 W19-W23 N20-N22
      5  263 1340 (N14E56 L=205)   410     50        2 ( 249/1312- 278/1371) -53.1 E57-E53 N11-N15
      6 1637 1252 (N13W39 L=300)  1320    210        2 (1563/1179-1727/1341) -30.3 W34-W46 N08-N17
      7  742 1302 (N12E18 L=242)   268     61        0 ( 715/1269- 766/1334) -36.9 E20-E16 N09-N13
      8 1203 1259 (N11W10 L=271)   190     91        0 (1156/1239-1245/1281)  21.8 W07-W12 N09-N12
      9 1823 1155 (N09W56 L=316)   399     60        1 (1804/1120-1836/1188)  90.0 W53-W57 N06-N10
     10 1867 1108 (N06W61 L=321)  1141    238        4 (1823/1054-1908/1165) -80.5 W55-W66 N03-N09
     11  154 1096 (S01E65 L=195)   668     84        0 ( 134/1054- 177/1141)  53.1 E68-E62 S03-N01
     12  387 1019 (S06E41 L=219)   948    191        0 ( 349/ 939- 461/1076)  58.0 E44-E35 S10-S02
     13  545  963 (S09E30 L=231)   306     93        0 ( 514/ 927- 589/ 997)  26.6 E32-E26 S11-S07
     14  736  798 (S19E17 L=243)  1105    206        2 ( 645/ 685- 868/ 878)  34.5 E24-E08 S25-S13
     15 1722  784 (S14W49 L=310)   434     82        1 (1679/ 755-1753/ 805)  -8.1 W45-W52 S16-S13
     16 1546  608 (S27W39 L=299)  1303    243        2 (1437/ 552-1620/ 710)  45.0 W31-W42 S30-S19
     17  599  543 (S35E30 L=231)   768    187        2 ( 544/ 481- 667/ 606)  33.7 E35-E23 S39-S31
     18 1768  565 (S27W60 L=321)   929    119        0 (1720/ 531-1808/ 593)  29.7 W55-W64 S29-S25
  

    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