Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  458 1665 (N58E54 L= 29)   856     88        0 ( 421/1636- 504/1688)  18.4 E61-E45 N56-N60
      2 1624 1424 (N19W50 L=133)   508    174        0 (1559/1371-1679/1455) -45.0 W44-W53 N14-N21
      3  643 1410 (N36E21 L= 62)   391    113        0 ( 581/1392- 699/1426)  -5.7 E26-E16 N34-N36
      4  283 1058 (N17E53 L= 30)   488     50        2 ( 266/1028- 305/1091) -90.0 E54-E50 N15-N19
      5 1883 1008 (S11W65 L=148)   711     69        1 (1847/ 985-1919/1033)   0.0 W59-W71 S12-S10
      6 1379  682 (S20W17 L=100)   345    114        1 (1339/ 633-1417/ 732)  50.2 W14-W20 S23-S17
      7 1056  570 (S22E05 L= 78)   956    272        0 ( 962/ 493-1145/ 673)  29.7 E12-W02 S26-S16
      8 1530  568 (S32W29 L=112)   396    128        0 (1480/ 517-1570/ 623)  40.6 W24-W31 S34-S28
      9  579  570 (S16E37 L= 46)   405     96        0 ( 550/ 526- 618/ 602)  39.8 E40-E33 S18-S13
     10 1573  509 (S37W33 L=116)   433     90        0 (1534/ 490-1612/ 543)  15.9 W29-W36 S38-S35
     11 1103  464 (S31E04 L= 79)   823    331        1 ( 962/ 395-1231/ 536)  14.0 E15-W05 S33-S27
     12  684  428 (S27E35 L= 48)   519     88        0 ( 662/ 392- 707/ 467)  71.6 E36-E32 S29-S23
  

    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