Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1253 1853 (N65W43 L=349)  1366    141        1 (1186/1831-1310/1885)   3.3 W31-W67 N62-N67
      2 1543 1574 (N37W47 L=354)   370     62        0 (1514/1558-1578/1589)   0.0 W44-W51 N36-N38
      3  935 1377 (N27E05 L=301)  1104    386        2 ( 874/1291-1043/1469) -54.0 E09-W02 N20-N33
      4 1541 1384 (N24W39 L=345)   282     65        0 (1517/1351-1560/1415) -76.0 W37-W40 N22-N26
      5 1455 1224 (N15W29 L=336)   292     67        0 (1430/1190-1484/1249)  36.9 W27-W31 N12-N16
      6 1167 1139 (N11W09 L=316)   198     38        1 (1140/1123-1192/1158) -18.4 W07-W10 N09-N12
      7  220 1119 (N10E61 L=245)   437     48        1 ( 207/1092- 240/1151) -63.4 E62-E58 N08-N12
      8 1746  857 (S10W51 L=357)   616      5        1 (1721/ 830-1772/ 882)  71.6 W48-W53 S11-S08
      9  774  814 (S08E16 L=290)   763    199        0 ( 697/ 774- 837/ 849)  -6.3 E21-E12 S10-S05
     10  322  654 (S19E55 L=251)   593     69        2 ( 290/ 627- 360/ 675)  18.4 E59-E50 S20-S17
     11 1531  596 (S26W35 L=342)   465    114        1 (1495/ 569-1576/ 636) -26.6 W31-W39 S27-S22
     12 1168  150 (S67W16 L=322)   721     40        1 (1150/ 121-1180/ 180) -53.1 W12-W25 S74-S62
  

    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