Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1255 1762 (N56W09 L=167)   576     57        1 (1227/1732-1284/1795)  58.0 W05-W12 N52-N60
      2 1524 1708 (N56W43 L=201)   452     38        1 (1508/1684-1546/1730)  29.1 W38-W47 N53-N58
      3  866 1528 (N30E19 L=140)   462     86        2 ( 837/1483- 889/1570) -63.4 E20-E16 N27-N33
      4  343 1455 (N18E58 L=100)  1914    274        2 ( 274/1351- 410/1557)  75.1 E62-E55 N10-N25
      5  263 1471 (N18E71 L= 88)  1343     98        5 ( 229/1434- 308/1508)   8.7 E77-E63 N15-N20
      6 1561 1407 (N32W34 L=193)  1128    283        1 (1441/1337-1660/1470) -14.7 W24-W43 N28-N34
      7  407 1404 (N16E49 L=109)   688    108        4 ( 373/1357- 443/1444)  68.2 E52-E46 N12-N18
      8  892  929 (S07E07 L=151)   311     40        1 ( 870/ 905- 918/ 953) -18.4 E08-E05 S08-S05
      9  198  899 (S18E62 L= 96)   900     80        1 ( 170/ 861- 226/ 937)  32.0 E66-E58 S20-S15
     10 1236  520 (S28W22 L=180)   736    295        0 (1089/ 472-1338/ 579)  38.2 W12-W28 S33-S22
     11 1569  521 (S24W48 L=207)   437     63        0 (1539/ 488-1591/ 548)  78.7 W46-W50 S26-S21
     12  443  441 (S47E46 L=113)   528     59        1 ( 413/ 422- 481/ 457)  15.3 E50-E39 S49-S45
  

    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