Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1633 1631 (N22W65 L=251)   427     49        0 (1600/1616-1657/1647)  -7.1 W60-W69 N20-N22
      2 1655 1461 (N13W54 L=240)   338     54        0 (1635/1437-1674/1489)  33.7 W52-W56 N11-N13
      3 1696 1428 (N10W56 L=242)   328     59        1 (1670/1408-1716/1451)   0.0 W52-W58 N08-N10
      4  977 1359 (N27W06 L=192)   477    167        1 ( 900/1336-1073/1378) -31.0 W01-W11 N23-N29
      5  541 1196 (N28E26 L=159)   261     82        1 ( 508/1174- 588/1210)   0.0 E28-E22 N27-N29
      6 1222 1184 (N11W15 L=201)   194     50        2 (1195/1168-1251/1204) -18.4 W13-W16 N09-N11
      7 1555 1147 (S00W34 L=220)   293    102        0 (1538/1098-1568/1195)  78.7 W32-W35 S03-N02
      8 1863 1092 (S14W59 L=245)   427     61        1 (1845/1066-1880/1118)  21.8 W56-W61 S15-S11
      9  205  899 (N17E61 L=124)  1117    253        2 ( 155/ 804- 290/ 998) -47.3 E66-E53 N10-N23
     10 1509  878 (S14W24 L=210)   244     74        2 (1462/ 853-1536/ 899) -11.3 W21-W26 S15-S12
     11 1803  844 (S26W48 L=234)   539     98        1 (1783/ 794-1821/ 893) -74.1 W46-W50 S29-S22
     12 1531  763 (S21W24 L=209)   394    141        0 (1477/ 722-1567/ 820)  23.2 W19-W26 S23-S18
     13  373  564 (S05E56 L=130)   381     83        2 ( 348/ 542- 402/ 597) -68.2 E57-E53 S07-S02
     14  920  419 (S26E24 L=162)  1689    328        2 ( 876/ 275- 983/ 570)  43.5 E33-E14 S35-S17
     15 1152  497 (S27E07 L=179)   474    122        0 (1113/ 453-1205/ 549)  20.6 E10-E02 S29-S24
     16 1692  496 (S47W37 L=223)   856     77        2 (1647/ 475-1726/ 525)  -5.7 W31-W43 S50-S44
     17  506  397 (S18E54 L=132)  1805    238        6 ( 437/ 318- 594/ 486) -79.7 E58-E47 S25-S12
  

    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