Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  740 1852 (N41E47 L= 21)  1115    199        2 ( 698/1807- 773/1883) -20.0 E52-E41 N36-N43
      2  787 1696 (N29E32 L= 36)   334     98        0 ( 770/1667- 807/1734) -45.0 E34-E30 N26-N31
      3  541 1439 (N08E39 L= 29)   640    220        3 ( 446/1395- 620/1470)   0.0 E46-E32 N06-N09
      4  529 1207 (S06E34 L= 35)  1322    360        3 ( 412/1048- 615/1364)  68.2 E39-E30 S16-N04
      5 1588 1191 (N15W31 L= 99)   480    176        0 (1564/1105-1621/1254) -69.4 W28-W34 N11-N19
      6 1162 1028 (S04W08 L= 76)   308    135        1 (1117/ 974-1221/1087)  66.8 W05-W09 S07-N00
      7  754 1023 (S13E16 L= 52)   187     57        2 ( 722/1007- 787/1038)   0.0 E17-E13 S13-S11
      8  296  963 (S24E48 L= 20)   450    109        1 ( 256/ 944- 330/ 987)  15.9 E51-E44 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