Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1522 1607 (N30W45 L= 60)   552    101        0 (1471/1568-1568/1635) -45.0 W41-W47 N27-N32
      2 1045 1471 (N26W06 L= 21)   481     99        0 ( 995/1446-1096/1493)   0.0 W02-W09 N24-N27
      3 1768 1460 (N18W62 L= 77)   456     49        0 (1739/1445-1801/1473)   6.3 W58-W67 N17-N18
      4  763 1340 (N21E13 L=  1)   259     63        0 ( 741/1313- 792/1374) -53.1 E14-E11 N18-N23
      5  339  607 (S18E55 L=320)   570     67        1 ( 306/ 587- 378/ 632)  -8.1 E58-E51 S19-S16
      6  389  521 (S24E54 L=320)   500    115        2 ( 360/ 466- 414/ 573)  28.6 E61-E50 S27-S21
      7 1691  418 (S48W60 L= 75)   571     60        0 (1663/ 390-1718/ 443) -14.7 W52-W71 S50-S45
      8 1629  389 (S50W52 L= 67)   464     49        0 (1598/ 368-1656/ 404) -14.0 W46-W59 S51-S48
      9 1521  349 (S52W40 L= 54)  1121    138        0 (1478/ 317-1566/ 395) -26.6 W32-W48 S55-S47
  

    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