Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  882 1861 (N43E47 L=184)   531     31        2 ( 863/1838- 900/1884) -45.0 E50-E42 N40-N44
      2 1245 1734 (N45E11 L=220)  1593    219        1 (1183/1662-1321/1783)  74.1 E13-E04 N37-N51
      3 1584 1393 (N33W25 L=255)   728    130        4 (1559/1332-1605/1455) -60.3 W22-W26 N29-N36
      4 1746 1371 (N38W41 L=272)  1434    228        2 (1667/1278-1804/1449) -21.8 W32-W48 N33-N41
      5 1348 1409 (N25W08 L=238)   247     61        0 (1317/1391-1385/1428)   0.0 W05-W10 N24-N26
      6  800  761 (S27E06 L=225)   386    141        1 ( 747/ 715- 862/ 797)  63.4 E08-E02 S30-S22
      7 1267  669 (S18W25 L=255)   436    212        0 (1184/ 644-1332/ 713)  32.0 W19-W28 S20-S13
      8 1002  461 (S38W18 L=249)   515     89        1 ( 967/ 432-1045/ 496)  74.1 W15-W20 S41-S34
      9 1150  336 (S41W36 L=267)   314     55        0 (1123/ 321-1176/ 358)  76.0 W34-W37 S42-S38
  

    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