Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1236 1828 (N58E14 L=219)  1007     87        1 (1205/1791-1258/1871) -36.0 E19-E08 N53-N64
      2  912 1722 (N37E30 L=203)   386    118        2 ( 889/1674- 935/1785) -32.5 E36-E25 N33-N41
      3 1451 1582 (N44W15 L=248)   513     75        0 (1418/1547-1477/1607) -36.9 W12-W17 N41-N45
      4 1434 1408 (N30W16 L=249)   271     74        1 (1419/1369-1449/1450) -51.3 W14-W18 N27-N32
      5 1800 1170 (N26W51 L=284)  2357    449        3 (1663/1004-1904/1339) -26.6 W36-W67 N18-N33
      6  353  976 (S21E44 L=189)   324     85        1 ( 326/ 948- 378/ 997)  -9.5 E46-E40 S23-S19
      7 1080  840 (S12W08 L=240)   195     66        0 (1056/ 808-1100/ 874)  90.0 W06-W08 S14-S09
      8 1266  494 (S28W31 L=264)  1808    392        3 (1094/ 431-1441/ 554)  41.6 W22-W40 S36-S19
      9 1494  443 (S24W51 L=283)   618    118        0 (1443/ 405-1538/ 481)  71.6 W48-W52 S27-S21
     10 1509  295 (S32W68 L=300)  1087    108        1 (1459/ 282-1567/ 319)   9.0 W60-W79 S34-S30
  

    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