Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1046 1763 (N48W06 L=282)  1059    131        3 (1000/1727-1101/1823) -59.0 W01-W10 N44-N54
      2 1510 1648 (N37W43 L=319)  1626    343        1 (1340/1608-1651/1678)  -5.9 W27-W58 N33-N39
      3 1274 1516 (N28W19 L=295)  3579    747        1 ( 943/1358-1553/1601) -19.3 E03-W37 N16-N34
      4  963 1385 (N20E03 L=273)   394    103        0 ( 925/1340-1011/1426) -56.3 E04-W00 N16-N22
      5 1834 1311 (N13W63 L=339)  2127    250        4 (1764/1249-1904/1393) -45.0 W53-W72 N09-N18
      6  588 1327 (N18E27 L=248)   588    141        2 ( 533/1281- 643/1364) -40.6 E31-E23 N14-N20
      7 1685 1290 (N12W46 L=322)   643    117        1 (1633/1251-1727/1323) -35.5 W42-W49 N09-N14
      8 1435 1255 (N10W27 L=303)   233     78        0 (1397/1234-1466/1277)  11.3 W24-W29 N08-N11
      9  727 1156 (N07E18 L=258)   361    100        1 ( 685/1118- 762/1194) -51.3 E20-E15 N04-N09
     10 1740  898 (S12W49 L=325)  2113    312        2 (1663/ 765-1814/1043)  63.4 W43-W55 S19-S03
     11  639  969 (S04E24 L=252)   196     45        0 ( 627/ 938- 652/ 998)  76.0 E24-E23 S06-S02
     12 1093  835 (S14W04 L=279)   289    110        0 (1063/ 780-1121/ 880)  56.3 W01-W05 S17-S11
     13 1382  771 (S19W22 L=298)   612    249        0 (1351/ 658-1436/ 868)  72.9 W20-W25 S26-S13
     14 1755  656 (S27W56 L=332)  4978    591        8 (1542/ 481-1915/ 836)  39.4 W40-W76 S38-S15
     15  292  655 (S21E58 L=218)   809    168        3 ( 265/ 576- 335/ 716)  42.0 E66-E51 S26-S17
     16 1078  602 (S29W02 L=278)   246     99        0 (1044/ 565-1115/ 647)  51.3 W00-W04 S31-S26
     17 1287  464 (S39W18 L=294)   407     79        1 (1258/ 437-1327/ 503) -45.0 W15-W21 S41-S36
     18  945  428 (S41E09 L=267)   492    104        0 ( 896/ 405- 984/ 457)  15.9 E12-E05 S42-S38
  

    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