Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1364 1509 (N31W32 L=321)  1275    366        1 (1197/1436-1511/1560) -28.8 W20-W40 N23-N36
      2  234 1378 (N35E67 L=222)  1245    106        2 ( 182/1352- 288/1397)  -2.5 E81-E58 N32-N35
      3  287 1314 (N30E56 L=233)   420     48        1 ( 261/1290- 307/1339) -26.6 E59-E52 N28-N32
      4  772 1315 (N26E13 L=276)   466    167        1 ( 694/1299- 853/1332) -10.3 E18-E07 N25-N27
      5 1412 1221 (N13W27 L=316)   473    111        1 (1388/1169-1436/1264)  51.3 W25-W29 N09-N15
      6 1070 1096 (N10W04 L=293)   183     50        1 (1043/1075-1100/1118) -33.7 W02-W05 N08-N10
      7 1875  862 (S18W64 L=353)  5083    397        7 (1810/ 665-1937/1002)  74.1 W58-W76 S30-S09
      8 1685  890 (S12W43 L=332)   838    278        3 (1614/ 799-1754/ 978)  31.6 W36-W49 S16-S07
      9  251  684 (S09E63 L=226)  7916    588       16 ( 171/ 481- 315/ 867)  47.4 E79-E52 S23-N02
     10  315  735 (S06E54 L=235)   342     52        0 ( 300/ 710- 335/ 763) -63.4 E54-E52 S07-S03
     11 1612  717 (S22W36 L=325)   363    101        0 (1572/ 699-1644/ 751)   9.5 W33-W39 S22-S19
     12 1073  718 (S13E01 L=288)   231     72        1 (1038/ 705-1116/ 738) -33.7 E02-W01 S14-S12
     13 1750  686 (S27W51 L=340)   399     83        1 (1723/ 656-1782/ 718)  31.0 W48-W54 S28-S24
     14 1019  569 (S23E06 L=283)  1589    154        5 ( 960/ 470-1133/ 705) -86.6 E09-W00 S30-S13
     15  799  483 (S26E23 L=266)   311     60        0 ( 770/ 459- 833/ 506)  11.3 E24-E19 S27-S24
     16  462  451 (S26E51 L=237)   376     70        2 ( 436/ 421- 486/ 475) -76.0 E53-E48 S27-S23
     17  800  418 (S31E25 L=264)   424     97        0 ( 758/ 390- 824/ 460)  33.7 E28-E21 S32-S27
     18 1239  354 (S42W07 L=296)   921    119        2 (1178/ 328-1296/ 394)   5.7 W01-W11 S44-S39
  

    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