Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  944 1837 (N42E39 L=134)   954    160        2 ( 870/1799-1015/1876)   0.0 E49-E31 N40-N43
      2  467 1648 (N16E59 L=114)  2307    365        1 ( 333/1608- 613/1690)   7.1 E78-E45 N11-N18
      3 1558 1575 (N45W20 L=193)   399    113        0 (1516/1531-1598/1611) -20.6 W15-W24 N43-N46
      4 1507 1534 (N40W16 L=189)   383     40        0 (1490/1505-1525/1563) -45.0 W14-W17 N37-N41
      5  831 1245 (N01E17 L=157)   688    234        1 ( 778/1146- 906/1328) -22.6 E21-E09 S03-N04
      6  231 1127 (S20E58 L=115)   480     79        1 ( 207/1098- 260/1166) -38.7 E61-E54 S21-S17
      7  322 1120 (S18E48 L=125)   388     79        0 ( 304/1084- 345/1155)  68.2 E49-E46 S20-S15
      8  162 1053 (S25E66 L=107)  2644    232        3 ( 126/ 970- 197/1143)  57.5 E73-E61 S30-S19
      9 1732 1103 (N20W42 L=215)   970    141        0 (1673/1057-1781/1142)   0.0 W36-W47 N18-N22
     10  904 1074 (S07E08 L=165)   298     75        0 ( 874/1039- 935/1107)  68.2 E08-E06 S09-S04
     11  526  823 (S31E26 L=147)   303    108        0 ( 477/ 787- 569/ 850)  40.6 E29-E22 S34-S28
     12 1093  691 (S22W14 L=187)   795    234        2 ( 997/ 637-1193/ 728)  54.0 W10-W18 S28-S17
     13 1366  637 (S17W32 L=205)  1269    279        1 (1230/ 605-1502/ 666)  36.3 W24-W40 S22-S11
     14  826  555 (S39W03 L=176)   746    141        1 ( 789/ 497- 875/ 606) -63.4 W00-W05 S42-S34
     15 1618  453 (S18W61 L=235)   484     60        0 (1595/ 428-1649/ 477)  76.0 W59-W63 S19-S15
     16 1542  442 (S21W55 L=229)   447     63        1 (1520/ 410-1565/ 474) -68.2 W53-W59 S23-S18
     17  834  434 (S47W09 L=182)   256     34        0 ( 813/ 411- 853/ 454) -21.8 W06-W11 S47-S45
     18  531  395 (S62E20 L=153)   774    135        0 ( 482/ 352- 577/ 448) -15.9 E27-E13 S65-S58
  

    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