Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1114 1722 (N45E14 L=311)   320     86        1 (1083/1699-1163/1746) -12.5 E17-E08 N43-N46
      2  380 1633 (N21E70 L=255)  1064     90        2 ( 359/1593- 413/1677) -40.6 E74-E65 N18-N24
      3  823 1590 (N28E28 L=297)   499    172        0 ( 743/1561- 907/1616)   4.8 E34-E21 N26-N28
      4  450 1455 (N11E49 L=276)   353     65        1 ( 429/1435- 487/1479)  26.6 E50-E45 N09-N12
      5 1171 1373 (N22W01 L=326)   323    175        1 (1119/1308-1235/1435) -29.1 E03-W06 N19-N24
      6 1658 1235 (N25W38 L=  2)   758    147        3 (1636/1160-1673/1306)  84.3 W36-W39 N19-N29
      7  196 1167 (S10E64 L=260)   519     37        1 ( 169/1146- 223/1180)   7.1 E68-E60 S11-S09
      8 1799  977 (N13W55 L= 19)   855     97        3 (1771/ 936-1836/1017)  56.3 W52-W58 N09-N15
      9  308  910 (S24E48 L=277)   635    142        2 ( 279/ 844- 330/ 974)  83.7 E50-E45 S28-S19
     10 1100  667 (S21W13 L=337)   844    300        0 (1062/ 550-1180/ 754)  85.9 W10-W17 S28-S14
     11  267  623 (S44E59 L=266)   638    101        0 ( 251/ 588- 293/ 676) -90.0 E62-E55 S46-S40
  

    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