Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  613 1439 (N37E25 L=127)   495     73        2 ( 580/1403- 645/1464)  31.0 E27-E22 N35-N39
      2  978 1372 (N28W02 L=154)   241     58        1 ( 948/1351-1014/1389) -14.0 W00-W04 N26-N29
      3  459  925 (N06E38 L=114)   770    210        2 ( 415/ 849- 504/1026)  51.3 E42-E33 N01-N11
      4  983  680 (S14E07 L=145)   255    126        0 ( 952/ 629-1003/ 743)  60.3 E09-E05 S17-S10
      5 1701  648 (S29W46 L=198)   720     59        1 (1671/ 617-1736/ 688)  45.0 W42-W49 S31-S26
      6  313  590 (S14E61 L= 92)   591    116        0 ( 296/ 543- 332/ 637)  40.6 E63-E55 S17-S11
      7  426  518 (S20E52 L=101)   407     64        2 ( 402/ 492- 461/ 547)   8.1 E54-E47 S21-S17
      8  641  370 (S32E40 L=112)   703    173        1 ( 597/ 314- 685/ 444)  31.0 E48-E33 S36-S26
      9  393  394 (S29E67 L= 86)   943     87        0 ( 373/ 358- 414/ 436)  19.2 E82-E59 S33-S25
  

    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