Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  414 1648 (N23E62 L= 89)  1004    225        0 ( 396/1597- 432/1697) -37.9 E67-E58 N19-N26
      2  452 1457 (N12E48 L=103)   296    113        0 ( 426/1420- 476/1494) -23.2 E51-E44 N09-N13
      3  336 1270 (S01E51 L=100)   509    141        1 ( 311/1216- 370/1327) -35.5 E54-E46 S04-N01
      4 1578 1203 (N21W31 L=183)   664    126        2 (1538/1139-1615/1264) -45.0 W27-W34 N16-N23
      5  620  901 (S18E22 L=129)   327    126        1 ( 580/ 868- 669/ 931)   0.0 E24-E18 S20-S15
      6  364  875 (S24E41 L=110)  1078    343        3 ( 278/ 798- 478/ 922)  32.0 E48-E32 S31-S20
      7 1131  831 (S11W10 L=162)   275     99        0 (1078/ 813-1180/ 852)  38.7 W07-W13 S13-S08
      8  905  715 (S23E01 L=151)   457    163        0 ( 833/ 686- 959/ 777)  58.0 E04-W01 S25-S17
      9 1416  447 (S27W42 L=193)  1039    215        0 (1330/ 414-1485/ 503)  45.0 W35-W45 S31-S22
     10  962  483 (S37W10 L=161)   454    146        0 ( 891/ 468-1041/ 501)  26.6 W04-W14 S38-S33
  

    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