Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1676 1597 (N21W69 L=186)   545     61        0 (1642/1577-1706/1610) -31.0 W65-W72 N19-N22
      2 1375 1551 (N28W39 L=156)   439     70        2 (1340/1524-1408/1576)   0.0 W35-W42 N27-N29
      3  251 1068 (N24E56 L= 62)   464     87        2 ( 232/1027- 274/1103) -78.7 E57-E52 N21-N26
      4  325 1007 (N20E48 L= 69)   677    193        2 ( 295/ 939- 389/1087) -60.9 E50-E42 N15-N24
      5 1749  916 (S18W45 L=163)   818    153        1 (1707/ 855-1801/ 981)  26.6 W40-W51 S21-S14
      6 1533  900 (S12W28 L=145)   516    134        0 (1471/ 885-1596/ 917) -18.4 W23-W32 S14-S10
      7  999  540 (S21E13 L=104)  2148    478        2 ( 836/ 397-1119/ 712)  25.7 E28-E01 S27-S12
      8 1235  496 (S30W00 L=118)   318     84        1 (1197/ 473-1267/ 528)  21.8 E02-W03 S31-S28
      9  724  396 (S25E37 L= 81)  1783    385        2 ( 657/ 327- 844/ 525) -86.4 E41-E28 S31-S15
     10 1622  355 (S59W40 L=158)   657     55        1 (1595/ 340-1654/ 374) -20.6 W33-W50 S61-S55
     11  655  269 (S33E52 L= 65)   619     99        2 ( 630/ 241- 677/ 311)  50.2 E57-E48 S35-S29
  

    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