Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1293 1712 (N42W42 L=252)  1426    273        3 (1165/1675-1414/1738)  -9.5 W29-W54 N38-N45
      2  746 1573 (N47E08 L=202)  2499    604        8 ( 550/1474-1038/1685)   0.0 E27-W19 N42-N49
      3 1539 1339 (N14W40 L=250)   666    221        1 (1462/1245-1590/1401) -71.6 W36-W41 N07-N19
      4  462 1344 (N36E35 L=175)   284     63        0 ( 430/1330- 497/1360) -15.9 E38-E31 N34-N37
      5  839 1116 (N16E09 L=201)   964    330        1 ( 758/ 986- 917/1258) -73.6 E12-E06 N07-N25
      6  548 1217 (N27E28 L=182)   265     83        1 ( 509/1204- 604/1233) -15.9 E30-E23 N25-N27
      7  187  970 (N16E63 L=147)   646     97        2 ( 165/ 926- 215/1017) -51.3 E66-E60 N13-N18
      8 1454  969 (S05W25 L=235)   326     82        0 (1424/ 941-1490/ 996) -59.0 W22-W26 S07-S02
      9 1523  850 (S14W27 L=237)   275     79        0 (1500/ 821-1545/ 883)  26.6 W25-W29 S15-S12
     10  667  748 (S03E27 L=183)   417    146        0 ( 651/ 684- 698/ 805)  45.0 E29-E23 S06-N00
     11  271  760 (N03E58 L=152)   604     74        1 ( 248/ 727- 297/ 802) -78.7 E59-E55 N00-N05
     12 1574  711 (S24W29 L=240)   713    147        1 (1525/ 668-1628/ 766) -45.0 W25-W33 S27-S19
     13  930  747 (S08E11 L=199)   227     80        0 ( 884/ 733- 964/ 761) -11.3 E13-E08 S08-S06
     14 1649  472 (S44W39 L=249)  2602    299        7 (1511/ 400-1772/ 544)   0.0 W23-W56 S47-S40
     15  348  438 (S19E68 L=142)  4599    276        7 ( 266/ 366- 422/ 515) -51.3 E78-E58 S24-S13
     16  567  399 (S24E47 L=163)  1213    214        3 ( 530/ 311- 596/ 492)  49.8 E52-E40 S30-S17
     17 1169  434 (S33E03 L=207)   249     82        0 (1137/ 406-1200/ 470)  26.6 E06-W00 S34-S30
     18 1432  346 (S48W16 L=226)   418     69        0 (1391/ 333-1471/ 358) -20.6 W12-W20 S50-S45
  

    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