Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  623 1802 (N30E55 L= 31)  1759    265        3 ( 593/1739- 655/1869) -50.7 E59-E49 N25-N36
      2 1613 1434 (N35W27 L=113)  3050    555        7 (1456/1226-1727/1593) -35.3 W12-W38 N24-N42
      3  416 1390 (S01E48 L= 38)  4060    703        8 ( 276/1112- 522/1566)  90.0 E61-E40 S19-N10
      4  283 1523 (N05E69 L= 17)  1165    329        4 ( 249/1495- 328/1562) -13.2 E78-E61 N03-N08
      5  900 1431 (N12E18 L= 68)   493    134        1 ( 838/1405- 987/1462)  40.6 E21-E14 N09-N16
      6 1297 1253 (N13W09 L= 95)  1936    413        2 (1140/1059-1397/1401) -30.6 E03-W19 N04-N17
      7  239 1252 (S12E59 L= 27)   776    140        1 ( 199/1192- 276/1315)  69.4 E61-E56 S15-S07
      8  548 1085 (S15E30 L= 57)   386    149        0 ( 502/1024- 617/1146)  60.9 E31-E26 S19-S10
      9  175 1040 (S25E60 L= 26)   928    229        1 ( 125/1012- 228/1066)  15.3 E65-E54 S27-S22
     10 1563  660 (S11W43 L=129)   288     90        0 (1550/ 615-1576/ 704) -51.3 W41-W45 S13-S08
     11 1560  447 (S22W54 L=140)   782    166        1 (1496/ 430-1625/ 469)  26.6 W48-W58 S23-S18
  

    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