Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1308 1645 (N35W31 L=143)  2503    601        8 (1049/1572-1577/1703)  -2.5 W09-W55 N31-N38
      2  307 1405 (N32E51 L= 61)   487     70        1 ( 271/1375- 332/1429) -24.0 E55-E46 N30-N34
      3  180 1265 (N25E63 L= 49)   879    154        1 ( 131/1199- 206/1325) -26.6 E77-E58 N20-N29
      4  410 1118 (N13E38 L= 74)   233     52        2 ( 382/1090- 431/1141)  45.0 E40-E36 N11-N13
      5  144 1035 (N11E66 L= 46)  3155    228       12 (  88/ 933- 196/1136) -81.3 E77-E58 N04-N17
      6 1636 1050 (S06W39 L=151)   394    115        0 (1621/ 993-1650/1112)  74.1 W37-W40 S09-S02
      7 1000  805 (S13E04 L=108)   319    143        0 ( 959/ 748-1028/ 850)  51.3 E07-E01 S15-S10
      8 1448  733 (S23W24 L=136)   947    283        1 (1381/ 616-1539/ 833)  45.0 W18-W30 S29-S17
      9 1734  813 (S22W47 L=159)   292     52        1 (1704/ 789-1762/ 833) -31.0 W43-W49 S23-S20
     10 1804  787 (S24W55 L=167)   599    105        2 (1769/ 757-1846/ 828)  23.2 W50-W59 S25-S21
     11  211  706 (S09E65 L= 47)   985    139        1 ( 187/ 645- 250/ 751)  71.6 E68-E59 S12-S06
     12  298  702 (S10E55 L= 57)   525    101        3 ( 260/ 669- 337/ 727) -38.7 E58-E51 S12-S08
     13  978  569 (S27E09 L=103)   452    168        1 ( 939/ 498-1007/ 658) -84.8 E10-E08 S32-S21
     14  444  615 (S17E45 L= 67)   321     49        0 ( 419/ 584- 466/ 641) -90.0 E46-E43 S19-S15
     15  517  583 (S20E41 L= 71)   277    113        0 ( 487/ 559- 556/ 610)   0.0 E42-E36 S21-S18
     16 1138  556 (S30W01 L=113)   580    154        0 (1075/ 511-1192/ 605)  32.0 E03-W05 S32-S27
     17  601  316 (S40E49 L= 63)  3165    425        3 ( 483/ 227- 761/ 477) -77.2 E56-E39 S49-S26
  

    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