Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  727 1685 (N25E39 L= 95)  1806    290        3 ( 691/1571- 778/1836) -53.1 E50-E31 N16-N36
      2  360 1579 (N11E63 L= 71)  1042    181        4 ( 313/1498- 400/1652) -57.5 E67-E59 N05-N16
      3 1114 1542 (N25E09 L=126)   952    255        1 (1069/1446-1143/1649) -50.2 E15-E05 N19-N31
      4  996 1532 (N21E15 L=119)   705    250        0 ( 970/1425-1031/1642) -68.2 E18-E11 N14-N29
      5 1764  979 (N11W48 L=182)   666    159        0 (1749/ 901-1777/1060) -73.3 W46-W50 N06-N16
      6  530  823 (S30E27 L=108)   577    153        0 ( 445/ 803- 599/ 839)   9.5 E33-E21 S32-S28
      7  624  671 (S38E17 L=117)   686    278        0 ( 501/ 646- 753/ 692)  24.2 E27-E07 S42-S32
      8  884  629 (S33W02 L=137)   394     79        1 ( 867/ 576- 903/ 670) -51.3 W00-W04 S36-S30
      9 1685  561 (S13W59 L=193)   597    115        0 (1631/ 536-1731/ 601)  54.5 W55-W61 S16-S09
     10  968  539 (S36W11 L=145)   707    142        1 ( 903/ 517-1036/ 561)  26.6 W06-W14 S38-S33
     11 1509  502 (S22W48 L=182)   369     94        1 (1468/ 477-1542/ 530)  68.2 W45-W50 S24-S19
  

    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