Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1397 1719 (N35W46 L=176)  1819    282        2 (1308/1657-1489/1762)  -7.1 W38-W54 N31-N36
      2  540 1604 (N50E24 L=105)   371     77        1 ( 519/1569- 556/1633) -59.0 E26-E22 N47-N52
      3  312 1370 (N38E45 L= 85)   524    115        1 ( 266/1351- 356/1388)   0.0 E49-E40 N37-N39
      4 1677 1191 (S02W44 L=174)  1417    306        2 (1629/1051-1741/1297)  34.5 W37-W53 S09-N02
      5 1393 1228 (N06W26 L=155)   327     81        1 (1352/1181-1427/1271) -71.6 W24-W27 N03-N09
      6  113  938 (N15E70 L= 59)   877    140        2 (  86/ 911- 142/ 964) -23.2 E74-E66 N12-N16
      7 1656  810 (S24W36 L=165)   424     81        0 (1613/ 790-1692/ 831)  -9.5 W32-W38 S25-S22
      8 1238  518 (S32W02 L=131)  1000    213        4 (1183/ 442-1266/ 611)  55.0 E03-W05 S36-S26
      9  795  432 (S28E29 L=100)   698    193        0 ( 765/ 368- 821/ 523)  51.3 E33-E25 S31-S21
     10 1490  480 (S42W20 L=150)   423     61        1 (1458/ 459-1515/ 502)   0.0 W17-W22 S44-S40
     11  366  448 (S18E62 L= 68)   600     88        0 ( 348/ 417- 386/ 476)  26.6 E64-E58 S19-S16
     12 1092  404 (S37E11 L=118)   283     62        0 (1065/ 383-1111/ 433)  36.9 E13-E09 S38-S34
  

    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