Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  458 1578 (N23E53 L= 51)   476     93        1 ( 433/1548- 488/1609) -33.7 E56-E49 N20-N24
      2 1669 1466 (N40W41 L=145)   524    115        1 (1629/1446-1715/1481)   7.1 W37-W46 N37-N41
      3 1445 1335 (N26W22 L=125)  1009    224        2 (1376/1229-1539/1404) -31.6 W16-W29 N20-N29
      4 1844 1058 (N16W58 L=162)  1150    194        0 (1810/ 986-1888/1119)  52.1 W55-W63 N11-N20
      5 1851  935 (N09W62 L=165)   543    163        3 (1818/ 901-1890/ 966)  35.5 W57-W66 N06-N11
      6 1682  781 (S03W48 L=152)   349    101        2 (1640/ 758-1713/ 809)  31.0 W45-W51 S04-S01
      7  318  742 (S31E46 L= 58)   633    137        0 ( 283/ 707- 373/ 787) -12.5 E49-E40 S32-S28
      8 1643  523 (S19W55 L=159)   930    151        2 (1599/ 489-1688/ 564)  71.6 W50-W59 S22-S16
  

    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