Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1284 1733 (N39W42 L=201)  1176    215        1 (1184/1679-1381/1770)  -2.7 W31-W52 N36-N40
      2 1449 1547 (N23W42 L=201)   331    104        0 (1413/1520-1496/1579) -68.2 W39-W43 N20-N25
      3  537 1347 (N33E25 L=134)   729    210        2 ( 448/1327- 651/1364) -14.0 E32-E16 N30-N35
      4 1137 1280 (N15W12 L=171)   256     76        0 (1089/1264-1174/1294) -31.0 W09-W14 N13-N16
      5 1815 1241 (S04W59 L=218)   397     64        0 (1787/1207-1836/1263) -45.0 W56-W60 S05-S01
      6  720 1242 (N22E13 L=146)   268     59        0 ( 690/1224- 755/1259)  11.3 E15-E10 N21-N23
      7 1755 1051 (S12W48 L=207)  1775    289        2 (1692/ 969-1817/1142) -63.4 W43-W55 S17-S05
      8  293 1014 (N18E49 L=110)  1786    345        2 ( 152/ 953- 417/1121)   3.8 E66-E36 N14-N21
      9 1851 1035 (S16W58 L=217)   614     47        1 (1831/1005-1874/1064)  26.6 W55-W60 S18-S14
     10 1588  830 (S21W31 L=190)   569    127        0 (1567/ 780-1604/ 894)  71.6 W28-W32 S24-S17
     11 1414  738 (S22W17 L=176)   804    225        0 (1347/ 646-1466/ 828)  38.7 W11-W21 S26-S17
     12 1203  655 (S22W02 L=161)   712    162        1 (1175/ 569-1230/ 732)  56.3 E01-W05 S27-S17
     13 1439  610 (S31W17 L=176)   503    131        0 (1384/ 581-1488/ 653)  14.0 W12-W20 S32-S29
     14  378  592 (S08E54 L=105)   320     33        0 ( 364/ 566- 393/ 616)  56.3 E55-E51 S09-S06
     15  849  446 (S27E26 L=133)   480    133        1 ( 814/ 388- 872/ 509)  37.9 E31-E22 S30-S23
     16 1644  457 (S50W37 L=196)  1106    120        2 (1578/ 435-1701/ 492)  -7.1 W28-W45 S52-S47
  

    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