Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  985 1646 (N30E17 L=290)   916    243        0 ( 926/1548-1048/1743) -38.7 E25-E10 N24-N36
      2  812 1676 (N29E30 L=278)   677    173        1 ( 766/1617- 880/1735)  84.3 E33-E26 N24-N34
      3  909 1638 (N28E22 L=286)   458    118        3 ( 878/1586- 946/1668)  90.0 E23-E19 N24-N31
      4  402 1596 (N17E58 L=249)  1214    130        3 ( 369/1547- 437/1651) -41.6 E64-E52 N13-N21
      5 1443 1623 (N41W16 L=323)   675    119        0 (1381/1586-1498/1651)   0.0 W10-W20 N39-N42
      6 1692 1509 (N40W40 L=348)   614     80        1 (1648/1478-1724/1533)  -7.1 W35-W44 N38-N41
      7  979 1338 (N10E09 L=299)   822    221        2 ( 939/1250-1040/1438)  90.0 E10-E07 N04-N17
      8 1158 1323 (N13W02 L=309)   327     86        0 (1138/1285-1178/1362) -68.2 W00-W03 N10-N15
      9 1364 1185 (N09W17 L=324)   584    217        0 (1298/1106-1409/1267)  74.7 W13-W18 N03-N14
     10 1027 1149 (N00E02 L=305)   453    154        0 ( 955/1127-1097/1190)  -7.1 E07-W01 S01-N01
     11 1587 1125 (N11W32 L=340)   250     44        0 (1564/1100-1608/1152) -33.7 W30-W33 N09-N12
     12 1841 1099 (N17W55 L=  2)   561    112        0 (1797/1073-1878/1120)  18.4 W50-W59 N14-N18
     13 1750 1049 (N11W46 L=354)   575     90        2 (1729/1002-1769/1096) -68.2 W44-W48 N08-N13
     14 1857  946 (N09W59 L=  7)  1128    161        6 (1812/ 873-1890/1020) -69.4 W54-W64 N04-N12
     15  842  961 (S14E09 L=298)   252    118        0 ( 796/ 928- 895/ 985)  45.0 E11-E06 S16-S11
     16 1508  871 (S05W32 L=340)   243     68        0 (1478/ 841-1528/ 892) -21.8 W29-W34 S06-S04
     17 1070  796 (S19W07 L=315)   529    185        1 (1012/ 724-1126/ 861)  68.2 W05-W09 S24-S14
     18 1132  459 (S38W22 L=330)   302     77        0 (1092/ 448-1180/ 474)  18.4 W19-W25 S39-S36
     19 1480  431 (S30W49 L=357)   291    133        0 (1422/ 411-1528/ 470)  50.2 W45-W50 S32-S26
  

    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