Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  942 1577 (N24E22 L= 23)   417    110        0 ( 902/1542-1000/1602)   0.0 E24-E16 N21-N24
      2 1543 1529 (N41W20 L= 64)   457     70        2 (1510/1506-1574/1563) -36.9 W16-W22 N39-N43
      3  384 1439 (N02E53 L=351)   324     60        2 ( 356/1420- 421/1457)   9.5 E56-E50 N00-N03
      4 1720 1316 (N33W38 L= 83)   405     88        0 (1682/1282-1743/1348) -18.4 W34-W41 N31-N34
      5 1775 1178 (N26W45 L= 90)   385     71        1 (1756/1143-1794/1208) -45.0 W43-W47 N24-N28
      6  601 1009 (S18E24 L= 20)   583    231        1 ( 506/ 964- 689/1083)  45.0 E30-E20 S22-S11
      7 1650 1018 (N12W38 L= 82)   246     68        0 (1635/ 979-1662/1050) -45.0 W35-W38 N10-N13
      8 1529 1019 (N08W29 L= 74)   307     96        1 (1499/ 994-1557/1043)  76.0 W27-W30 N06-N10
      9 1864  962 (N17W60 L=104)   614     77        2 (1838/ 921-1886/ 986)  90.0 W56-W62 N15-N19
     10  417  886 (S29E36 L=  9)   312     84        1 ( 382/ 856- 455/ 907)  38.7 E38-E32 S31-S27
     11 1482  299 (S31W65 L=109)   942     90        3 (1449/ 270-1524/ 339)   0.0 W57-W73 S32-S29
     12  949  303 (S51W25 L= 70)   315     59        0 ( 933/ 279- 972/ 331) -18.4 W22-W28 S52-S49
  

    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