Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1052 1520 (N36W15 L=296)  1351    336        1 ( 921/1488-1212/1544) -19.3 W05-W25 N31-N39
      2 1435 1485 (N24W39 L=320)   988    287        1 (1325/1412-1537/1532) -50.2 W32-W43 N16-N28
      3 1854 1186 (S05W63 L=344)  2141    197       10 (1796/1130-1922/1256) -37.6 W56-W72 S10-N00
      4  798 1113 (N17E12 L=269)   617    125        3 ( 773/1036- 832/1163)  82.9 E13-E09 N12-N20
      5  333 1034 (N20E47 L=234)   748    168        0 ( 296/ 982- 396/1095) -45.0 E50-E41 N16-N24
      6  324  861 (N10E50 L=231)   294     75        0 ( 308/ 826- 342/ 905) -78.7 E51-E48 N07-N12
      7 1761  784 (S26W46 L=327)  1333    250        2 (1680/ 680-1815/ 880)  26.6 W37-W52 S31-S21
      8  358  566 (S08E57 L=224)  6913     15       19 ( 247/ 386- 509/ 757)  30.3 E74-E38 S21-N02
  

    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