Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  682 1324 (N01E27 L=359)   512    118        0 ( 632/1291- 737/1356)  -7.1 E31-E23 N00-N02
      2 1595 1113 (N14W31 L= 58)   638    236        0 (1548/1009-1622/1205) -38.7 W25-W35 N08-N18
      3 1492 1158 (N13W23 L= 50)   282     73        1 (1475/1113-1505/1191) -68.2 W22-W24 N10-N15
      4 1727 1006 (N12W44 L= 70)   317    105        0 (1704/ 960-1744/1044) -36.9 W41-W45 N10-N14
      5 1792  779 (N03W57 L= 84)   475     84        2 (1758/ 747-1824/ 814)  68.2 W53-W59 N00-N05
      6  513  771 (S35E27 L=360)   799    164        3 ( 441/ 717- 576/ 807)  -4.8 E33-E21 S36-S32
      7 1651  465 (S18W63 L= 89)   401     38        1 (1630/ 446-1668/ 488)  63.4 W60-W65 S18-S16
      8  477  384 (S65E33 L=354)  1584    227        4 ( 379/ 344- 572/ 419)   4.0 E60-E16 S68-S61
  

    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