Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  984 1772 (N38E30 L= 64)   523     73        1 ( 947/1743-1026/1798)  -7.1 E34-E26 N36-N39
      2 1174 1537 (N28E07 L= 88)   846    200        1 (1101/1459-1216/1640) -42.3 E14-E02 N23-N33
      3  917 1414 (N13E17 L= 77)   279     64        0 ( 888/1396- 950/1435)  33.7 E18-E14 N11-N14
      4 1506 1300 (N24W21 L=115)   419     85        0 (1462/1267-1538/1326)   9.5 W17-W23 N22-N25
      5 1151 1290 (N12E00 L= 94)   474    165        0 (1065/1252-1215/1324)   0.0 E05-W04 N11-N13
      6 1469 1149 (N14W22 L=116)   260     46        1 (1452/1113-1485/1177) -76.0 W20-W22 N11-N15
      7 1695 1060 (N16W40 L=134)   849    158        1 (1639/1007-1758/1110)  -4.8 W34-W46 N14-N17
      8 1588 1004 (N09W33 L=128)   244     76        0 (1559/ 974-1615/1029)  -9.5 W30-W36 N08-N10
      9 1098  577 (S28W19 L=113)   336     85        0 (1058/ 562-1149/ 594)  38.7 W16-W21 S30-S26
     10  549  523 (S50E19 L= 75)   385     68        0 ( 514/ 511- 593/ 539)  18.4 E23-E14 S51-S48
     11 1528  391 (S24W59 L=153)  2704    469        4 (1359/ 286-1657/ 493)  62.2 W51-W64 S34-S15
     12 1115  417 (S37W28 L=122)   399     97        0 (1070/ 377-1154/ 456)  81.9 W26-W29 S41-S33
  

    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