Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  649 1732 (N40E48 L=166)   519     54        1 ( 624/1710- 680/1758)  -7.1 E51-E43 N38-N42
      2 1508 1674 (N53W35 L=249)   463     68        0 (1473/1659-1548/1687)  14.0 W29-W41 N51-N55
      3 1359 1600 (N43W16 L=231)   278     69        0 (1327/1577-1387/1618) -21.8 W13-W18 N41-N44
      4 1648 1401 (N33W42 L=256)   359     88        0 (1605/1386-1698/1418)   0.0 W37-W47 N32-N34
      5  450  665 (S31E37 L=178)   385    135        0 ( 396/ 632- 506/ 696)  28.6 E42-E31 S34-S27
      6 1709  648 (S13W56 L=270)   370     31        1 (1693/ 624-1726/ 670) -36.9 W53-W58 S14-S11
      7 1530  460 (S28W49 L=263)   576    109        0 (1475/ 443-1587/ 477)  11.3 W43-W54 S29-S26
      8 1166  362 (S41W24 L=239)   330     66        2 (1127/ 345-1197/ 383)   0.0 W21-W27 S41-S39
      9  747  211 (S67E13 L=202)   690     42        1 ( 721/ 186- 779/ 238) -50.2 E17-E07 S69-S62
  

    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