Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1425 1802 (N43W58 L= 77)   616    122        1 (1402/1770-1444/1825)  31.0 W55-W61 N41-N44
      2 1419 1495 (N25W38 L= 56)   362     99        0 (1393/1455-1452/1528)  15.9 W34-W41 N23-N26
      3 1043 1450 (N32W11 L= 30)   397    108        0 ( 994/1427-1106/1465) -14.0 W07-W15 N30-N32
      4  422 1376 (N40E38 L=341)  1164    175        2 ( 365/1328- 485/1420)  16.7 E43-E32 N36-N41
      5  636 1141 (N21E22 L=357)   280     64        0 ( 601/1121- 674/1156) -11.3 E24-E19 N19-N22
      6  905 1062 (N12E06 L= 12)   585    110        1 ( 841/1037- 962/1090) -33.7 E09-E03 N09-N14
      7  261  944 (N14E55 L=324)   564    115        1 ( 244/ 889- 280/1006)  82.9 E56-E52 N10-N18
      8  464  790 (N03E40 L=338)   308     51        0 ( 433/ 769- 493/ 815) -45.0 E42-E38 N01-N04
      9  675  490 (S19E34 L=344)   412    127        1 ( 658/ 443- 699/ 538)  45.0 E36-E30 S22-S16
     10 1631  454 (S46W36 L= 54)   577    103        0 (1593/ 435-1665/ 470) -14.0 W31-W39 S48-S43
     11  610  262 (S35E52 L=326)   600     87        1 ( 591/ 237- 641/ 284)  18.4 E56-E49 S37-S33
     12 1076  180 (S54E23 L=355)   438     58        0 (1040/ 172-1109/ 195)  -8.1 E26-E19 S55-S52
  

    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