Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1619 1666 (N32W62 L= 29)   664    210        0 (1582/1639-1641/1704)   7.1 W57-W65 N31-N33
      2 1418 1543 (N30W39 L=  5)  1158    203        2 (1340/1499-1504/1591)   0.0 W32-W47 N28-N32
      3  495 1485 (N44E34 L=292)   402     77        3 ( 457/1464- 531/1505)   0.0 E37-E30 N42-N45
      4 1573  949 (S08W33 L=  0)  1348    261        5 (1521/ 833-1607/1076)  58.4 W28-W37 S14-N00
      5 1773 1045 (S07W51 L= 18)   434     94        1 (1746/1016-1800/1073) -45.0 W49-W54 S08-S04
      6  234  920 (N10E58 L=268)  1154    220        2 ( 180/ 866- 291/ 990)  23.2 E65-E51 N06-N13
      7 1693  930 (S12W42 L=  9)   649    127        1 (1670/ 866-1721/ 987)  60.3 W39-W45 S15-S08
      8 1853  897 (S18W58 L= 25)   551     87        1 (1839/ 854-1866/ 937)  68.2 W56-W59 S20-S15
      9 1319  437 (S37W10 L=337)  1359    296        1 (1227/ 355-1390/ 564)  40.6 W02-W16 S42-S28
  

    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