Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  565 1678 (N20E51 L=132)  1704    316        6 ( 507/1616- 636/1744) -11.3 E60-E45 N15-N23
      2 1821  935 (N14W54 L=237)   403     64        2 (1799/ 909-1840/ 962)  63.4 W51-W55 N11-N15
      3  851  905 (S17E06 L=177)   605    162        1 ( 766/ 854- 925/ 952)  55.0 E10-E02 S22-S12
      4 1584  734 (S05W42 L=225)  1499    351        4 (1546/ 601-1619/ 913) -60.6 W36-W46 S12-N04
      5  213  880 (S35E54 L=129)   503     79        2 ( 189/ 851- 232/ 912)  76.0 E56-E52 S36-S32
      6 1150  672 (S22W18 L=201)   349     98        0 (1102/ 650-1199/ 702)   8.1 W14-W21 S22-S20
      7 1511  618 (S13W42 L=225)   264    105        0 (1465/ 586-1532/ 639)  -8.1 W37-W44 S14-S12
      8  653  575 (S43E11 L=172)  1319    356        2 ( 493/ 552- 816/ 604)  19.1 E25-W01 S48-S38
      9 1490  384 (S26W54 L=237)  1368    310        4 (1360/ 355-1602/ 416)  30.5 W45-W62 S31-S20
  

    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