Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  570 1715 (N28E51 L=287)  3109    444        5 ( 482/1612- 640/1811) -79.4 E59-E43 N20-N36
      2  558 1493 (N13E39 L=298)  1213    265        2 ( 477/1415- 638/1601)  74.1 E42-E37 N07-N21
      3  697 1465 (N13E29 L=309)   388    117        1 ( 671/1414- 727/1508) -80.5 E30-E27 N10-N16
      4 1517 1432 (N28W24 L=  1)   268     70        1 (1496/1397-1542/1463) -53.1 W21-W25 N26-N30
      5  864 1314 (N07E15 L=323)   248     87        0 ( 838/1278- 898/1355)  68.2 E15-E13 N04-N09
      6 1819 1229 (N24W52 L= 30)   508     76        0 (1795/1188-1841/1269) -68.2 W49-W54 N21-N26
      7 1515 1231 (N15W26 L=  3)   389    104        1 (1469/1192-1544/1259) -18.4 W22-W28 N13-N16
      8 1334 1215 (N10W14 L=352)   281     66        2 (1312/1181-1356/1255) -53.1 W12-W15 N08-N12
      9  184 1065 (S17E63 L=275)  1909    232        3 ( 147/ 965- 238/1148)  63.4 E68-E57 S23-S11
     10  612 1093 (S10E26 L=311)   201     67        1 ( 574/1072- 635/1115)  26.6 E28-E24 S11-S08
     11 1767 1042 (N11W48 L= 26)   616    161        1 (1742/ 975-1816/1105)  66.0 W45-W52 N06-N15
     12 1544 1057 (N06W30 L=  8)   222     55        1 (1514/1037-1569/1084)  14.0 W28-W32 N04-N07
     13  548  820 (S27E27 L=310)  1307    537        3 ( 371/ 689- 655/1016)  50.0 E43-E21 S38-S13
     14 1779  678 (S09W61 L= 39)   449     72        1 (1737/ 657-1802/ 705)  38.7 W57-W63 S10-S06
     15 1575  517 (S23W51 L= 29)  1538    299        2 (1521/ 390-1624/ 651) -62.1 W47-W56 S31-S14
     16 1744  515 (S18W69 L= 46)  1503    183        2 (1673/ 453-1814/ 567)  42.3 W64-W76 S23-S13
  

    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