Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1640 1353 (N33W32 L=124)   552     71        1 (1618/1318-1663/1393) -45.0 W30-W34 N30-N35
      2 1815 1097 (N23W52 L=145)  2617    331        3 (1726/ 964-1913/1178)  40.6 W45-W67 N12-N30
      3  932 1104 (S02E07 L= 85)   249     59        0 ( 900/1079- 953/1130) -14.0 E09-E05 S03-S01
      4  412  765 (S35E36 L= 57)  2345    330        5 ( 324/ 658- 495/ 900)  52.7 E46-E29 S45-S24
      5 1767  647 (S04W63 L=155)   567     38        0 (1746/ 618-1783/ 675) -45.0 W59-W65 S05-S02
      6 1090  545 (S30W19 L=111)  1206    400        4 ( 949/ 458-1239/ 618)  56.3 W12-W25 S39-S21
      7 1376  579 (S19W35 L=127)   371    119        0 (1330/ 561-1425/ 616)  56.3 W31-W37 S21-S15
      8 1471  552 (S18W43 L=135)   518    111        0 (1426/ 517-1515/ 578)  63.4 W40-W45 S21-S15
      9 1647  499 (S15W60 L=152)   431     60        0 (1615/ 479-1678/ 523)  59.0 W58-W62 S17-S12
     10  696  479 (S47E06 L= 86)   251     81        0 ( 676/ 447- 712/ 518) -68.2 E07-E04 S48-S43
  

    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