Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1761 1468 (N18W65 L=243)  1017    236        1 (1713/1448-1812/1489) -13.0 W59-W72 N15-N19
      2 1356 1329 (N19W27 L=205)  1038    316        2 (1231/1261-1517/1395) -36.9 W19-W35 N11-N24
      3 1844 1201 (N01W63 L=242)   714    144        1 (1809/1174-1881/1227) -20.6 W59-W68 S01-N02
      4 1301 1096 (N06W18 L=196)   402    138        0 (1275/1034-1329/1173) -90.0 W16-W18 N02-N11
      5  836  955 (N05E12 L=166)   380    107        2 ( 790/ 901- 862/ 987) -71.6 E14-E10 N01-N07
      6 1737  836 (S18W45 L=224)   925    209        1 (1711/ 746-1757/ 953)  81.3 W42-W47 S24-S10
      7  692  885 (N03E22 L=156)   327     78        0 ( 673/ 839- 711/ 923) -90.0 E23-E21 N00-N05
      8  475  770 (S01E40 L=139)   394     95        1 ( 442/ 749- 519/ 791) -11.3 E41-E35 S02-N00
      9  880  499 (S24E19 L=160)   697    169        0 ( 802/ 474- 955/ 540)  -5.2 E23-E12 S25-S21
  

    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