Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1264 1626 (N35E02 L= 33)   278     27        0 (1245/1604-1286/1649) -33.7 E03-W00 N34-N36
      2 1372 1320 (N18W13 L= 47)   227     69        0 (1342/1285-1399/1351) -21.8 W10-W15 N17-N19
      3  200 1140 (S17E60 L=335)  1180    219        1 ( 170/1032- 228/1253)  90.0 E61-E58 S23-S10
      4  433 1144 (S13E39 L=355)  1039    243        2 ( 396/1052- 472/1239) -80.5 E42-E35 S19-S07
      5  499 1161 (S11E35 L=360)   297     46        0 ( 477/1131- 526/1190)  63.4 E36-E33 S13-S08
      6 1885 1003 (N16W59 L= 94)  2568    420        7 (1811/ 917-1934/1076) -10.0 W49-W67 N12-N20
      7  874  934 (S16E06 L= 28)   490    167        0 ( 816/ 855- 908/1006) -33.7 E11-E02 S19-S12
      8 1856  949 (N12W57 L= 92)   392     67        0 (1839/ 917-1872/ 979)  90.0 W55-W59 N10-N14
      9 1523  874 (S02W33 L= 67)  1004    224        1 (1430/ 826-1613/ 951)  -3.6 W24-W40 S04-N00
     10 1040  814 (S19W06 L= 41)   229     30        1 (1025/ 785-1055/ 842) -76.0 W05-W07 S20-S16
     11  570  782 (S32E23 L= 12)  1242    368        1 ( 409/ 751- 738/ 813)  11.3 E35-E10 S35-S28
     12 1868  768 (N03W67 L=101)   501     61        0 (1850/ 737-1883/ 801)  90.0 W65-W69 N01-N05
     13 1697  713 (S06W51 L= 86)   363    121        0 (1653/ 676-1716/ 745)  -8.1 W46-W54 S07-S04
     14 1588  593 (S15W48 L= 82)  1121    225        1 (1485/ 555-1680/ 640)  31.6 W40-W53 S18-S09
     15 1740  508 (S14W67 L=101)  3256    338        7 (1679/ 438-1802/ 592) -73.3 W60-W74 S19-S09
     16 1645  487 (S19W58 L= 93)   396     78        1 (1625/ 459-1671/ 513) -21.8 W55-W61 S19-S17
     17 1274  475 (S31W33 L= 67)   844    135        0 (1203/ 452-1356/ 506)  48.8 W28-W36 S34-S26
     18 1536  441 (S25W52 L= 87)   647    153        2 (1502/ 389-1561/ 501) -69.4 W48-W55 S28-S20
  

    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