Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  788 1416 (N16E19 L= 85)   207     70        0 ( 766/1383- 811/1447) -45.0 E20-E16 N14-N18
      2 1776 1172 (N13W50 L=154)   736    219        3 (1726/1095-1824/1221)  56.3 W45-W54 N07-N16
      3 1004  986 (S08E01 L=103)   825    274        1 ( 908/ 948-1120/1083)  -9.5 E07-W05 S09-S02
      4  567  991 (S11E28 L= 76)   346     80        0 ( 549/ 948- 593/1035) -59.0 E29-E26 S13-S08
      5  247 1000 (S12E54 L= 50)   341     83        0 ( 220/ 964- 266/1025) -53.1 E56-E51 S14-S10
      6  231  962 (S14E55 L= 49)   397    122        0 ( 201/ 934- 270/ 993) -23.2 E58-E51 S15-S12
      7  528  913 (S16E31 L= 74)   205     78        0 ( 500/ 888- 566/ 929)  31.0 E32-E27 S17-S14
      8 1562  863 (S09W36 L=140)   218     67        0 (1537/ 828-1582/ 896) -36.9 W33-W37 S10-S07
      9  921  728 (S24E04 L=101)   326    117        1 ( 904/ 665- 936/ 784)  90.0 E04-E02 S28-S20
     10 1119  672 (S26W10 L=114)   190     42        1 (1100/ 649-1144/ 690)  33.7 W08-W11 S26-S24
     11 1652  498 (S29W54 L=158)   412     90        0 (1607/ 466-1694/ 530)  45.0 W50-W56 S31-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