Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  692 1683 (N36E41 L= 65)   393     54        1 ( 674/1661- 717/1714) -15.9 E44-E37 N34-N37
      2  780 1459 (N22E24 L= 82)   326     95        0 ( 740/1422- 806/1495) -26.6 E27-E21 N19-N23
      3 1363 1267 (N20W18 L=123)   265     86        0 (1327/1243-1401/1294) -11.3 W14-W20 N18-N20
      4 1851 1212 (N25W62 L=168)  1803    152        3 (1779/1179-1927/1249)   0.0 W52-W79 N23-N26
      5 1179  651 (S20W16 L=122)   282    143        0 (1167/ 582-1200/ 714) -76.0 W15-W18 S24-S16
      6  928  513 (S34W03 L=108)   593    125        1 ( 903/ 458- 967/ 572) -54.5 W00-W06 S37-S30
      7 1272  477 (S30W28 L=133)   309     72        0 (1233/ 456-1302/ 508)  68.2 W25-W29 S32-S27
      8 1290  312 (S42W39 L=145)   686     87        0 (1237/ 288-1330/ 330)   0.0 W33-W44 S43-S40
  

    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