Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1141 1609 (N35E09 L= 76)  1155    207        3 (1099/1516-1180/1694) -45.0 E15-E03 N29-N40
      2  283 1506 (N08E70 L= 15)   720     76        2 ( 265/1478- 311/1538) -26.6 E75-E64 N06-N10
      3  520 1465 (N10E44 L= 41)   353    116        0 ( 503/1413- 543/1518) -71.6 E45-E41 N07-N13
      4 1449 1233 (N19W20 L=105)   890    205        4 (1395/1150-1499/1328)  72.9 W18-W22 N13-N26
      5 1887 1265 (N35W66 L=151)   617     70        0 (1876/1225-1898/1300)  40.6 W63-W71 N32-N38
      6 1697 1098 (N18W41 L=126)  1040    221        2 (1666/1006-1733/1200) -48.0 W37-W46 N13-N23
      7 1113 1112 (N04W03 L= 88)   194     56        0 (1081/1098-1143/1128)  45.0 W01-W04 N02-N05
      8 1816  920 (N11W57 L=142)   726     95        1 (1786/ 887-1850/ 959)  63.4 W54-W61 N08-N14
      9  367  771 (S34E41 L= 44)   738    183        2 ( 288/ 735- 437/ 821)  -7.1 E49-E33 S35-S31
     10  480  785 (S30E30 L= 55)   286     66        2 ( 451/ 765- 504/ 807)  45.0 E32-E28 S32-S28
     11 1480  620 (S15W39 L=124)   311    143        2 (1433/ 575-1536/ 660)  56.3 W37-W42 S17-S11
  

    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