Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1294 1730 (N54W17 L= 84)   374     51        0 (1263/1717-1328/1743)   0.0 W13-W20 N52-N55
      2  585 1475 (N25E37 L= 30)  1350    323        3 ( 490/1393- 635/1608) -82.4 E42-E32 N19-N34
      3 1474 1406 (N30W29 L= 95)  2306    598        2 (1223/1235-1605/1530) -27.3 W09-W38 N19-N37
      4 1656 1231 (N20W43 L=109)   756    166        2 (1620/1174-1690/1296) -74.1 W39-W45 N16-N23
      5 1334 1201 (N15W18 L= 85)   351    146        0 (1302/1145-1380/1269) -49.4 W15-W21 N12-N19
      6 1520 1181 (N15W31 L= 98)   323    120        0 (1503/1126-1531/1243) -81.9 W30-W32 N12-N19
      7  423  781 (S20E39 L= 27)   658    149        1 ( 401/ 715- 454/ 843)  76.0 E41-E36 S23-S15
      8 1434  650 (S18W31 L= 98)   715    200        2 (1353/ 585-1507/ 711)  45.0 W26-W35 S22-S13
      9 1708  652 (S16W55 L=121)   568     52        1 (1685/ 624-1736/ 679)  63.4 W52-W57 S17-S13
     10 1720  577 (S21W60 L=126)  1343    193        3 (1660/ 514-1780/ 633)  56.3 W53-W68 S25-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