Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  724 1698 (N48E30 L=342)   368     90        0 ( 708/1667- 740/1731) -39.8 E33-E27 N45-N50
      2  679 1619 (N41E31 L=342)   451     80        1 ( 647/1593- 704/1654)  63.4 E32-E28 N38-N43
      3 1065 1478 (N32W02 L= 14)   462    101        1 (1020/1443-1112/1511) -33.7 E01-W05 N29-N34
      4 1328 1418 (N29W20 L= 33)   284     73        1 (1305/1385-1353/1443)  45.0 W18-W22 N26-N30
      5 1577 1325 (N23W38 L= 51)  1920    410        5 (1415/1238-1718/1415) -25.6 W26-W49 N16-N28
      6  627 1328 (N21E28 L=345)   572    155        1 ( 605/1242- 660/1389) -66.0 E29-E24 N15-N24
      7 1163 1273 (N19W08 L= 21)   974    256        3 (1105/1127-1216/1374) -69.4 W04-W11 N09-N25
      8 1320 1264 (N18W19 L= 32)   744    162        2 (1258/1189-1393/1342) -51.3 W14-W23 N13-N23
      9  970 1062 (N05E03 L=  9)  1276    360        1 ( 884/ 970-1120/1200) -38.2 E09-W05 N00-N13
     10  821  921 (S04E12 L=  1)   336    118        1 ( 772/ 903- 884/ 939)   8.1 E15-E08 S05-S02
  

    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