Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  517 1653 (N43E51 L= 41)   563     52        0 ( 490/1624- 541/1678) -15.3 E57-E46 N40-N44
      2 1397 1572 (N40W30 L=122)   387     62        1 (1364/1545-1424/1594) -36.9 W27-W31 N37-N41
      3  290 1515 (N31E71 L= 21)   855    131        0 ( 273/1460- 306/1573) -21.3 E81-E63 N27-N34
      4 1669 1448 (N31W52 L=144)   656    117        0 (1628/1405-1705/1488) -56.3 W48-W55 N27-N33
      5 1198 1369 (N25W11 L=103)   346    100        0 (1160/1331-1245/1402)  33.7 W08-W14 N22-N27
      6 1449 1263 (N19W28 L=121)   561    216        0 (1369/1188-1511/1327) -45.0 W22-W32 N13-N22
      7 1850 1260 (N18W67 L=160)  1349     98        3 (1823/1207-1873/1314)  74.1 W64-W70 N14-N21
      8 1814 1200 (N14W60 L=153)   821    121        4 (1773/1160-1848/1244) -63.4 W55-W64 N11-N17
      9 1052  873 (S06W02 L= 95)   465    233        0 ( 987/ 822-1090/ 950)  69.4 E01-W04 S09-S01
     10 1748  855 (S07W52 L=145)   511     44        1 (1724/ 831-1774/ 877)  36.9 W49-W54 S08-S05
     11  885  749 (S14E08 L= 84)  1128    268        3 ( 847/ 648- 923/ 873)  82.4 E10-E05 S21-S06
     12  440  786 (S14E40 L= 53)   345     59        2 ( 413/ 767- 476/ 805)  11.3 E41-E36 S15-S12
     13 1185  670 (S19W11 L=104)   340    120        0 (1141/ 632-1218/ 712) -45.0 W08-W13 S21-S16
     14  269  686 (S22E59 L= 33)   459     67        0 ( 243/ 665- 293/ 710) -33.7 E62-E55 S22-S20
     15 1454  523 (S29W33 L=126)  1523    245        3 (1411/ 431-1497/ 634) -59.0 W30-W40 S35-S20
     16  746  599 (S25E18 L= 74)   638    173        3 ( 677/ 572- 820/ 633)  -5.7 E23-E12 S26-S22
     17 1669  494 (S31W57 L=149)   557     80        0 (1626/ 481-1709/ 507)   0.0 W50-W62 S32-S30
  

    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