Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1437 1461 (N20W38 L=249)   619    136        1 (1372/1446-1502/1481) -24.0 W33-W42 N17-N21
      2 1411 1325 (N13W31 L=242)   239     42        2 (1387/1307-1437/1345) -14.0 W28-W32 N12-N14
      3  660 1195 (N25E18 L=193)   233     90        1 ( 617/1178- 708/1213) -33.7 E20-E14 N23-N27
      4  225 1105 (N29E57 L=154)  1189    205        3 ( 139/1052- 294/1174) -22.2 E71-E49 N24-N33
      5  600 1010 (N16E25 L=185)   227     87        0 ( 568/ 975- 629/1039) -59.0 E27-E24 N13-N18
      6  852  998 (N10E10 L=201)   238     69        0 ( 820/ 969- 880/1027) -68.2 E11-E08 N07-N12
      7  454  916 (N14E38 L=173)   351    177        2 ( 430/ 841- 472/ 987)  76.0 E40-E35 N09-N17
      8 1549  689 (S27W24 L=235)   518    111        0 (1514/ 641-1577/ 735) -63.4 W21-W25 S30-S22
      9 1840  690 (S38W55 L=266)   947    110        2 (1795/ 659-1893/ 720)  -3.6 W47-W65 S39-S36
     10 1072  637 (S16E07 L=204)   847    296        1 (1015/ 544-1167/ 705)  28.6 E11-W00 S21-S12
     11 1178  593 (S21E02 L=209)   405    117        1 (1132/ 554-1209/ 645)  31.0 E05-W00 S23-S18
     12  248  573 (S03E70 L=141)   638     48        1 ( 224/ 549- 274/ 598) -56.3 E73-E65 S05-S01
     13  994  527 (S21E15 L=196)   455     83        1 ( 951/ 507-1041/ 555) -38.7 E17-E12 S23-S19
     14  623  423 (S19E43 L=168)   300     36        0 ( 605/ 403- 645/ 442)  14.0 E44-E40 S20-S17
     15 1498  394 (S47W16 L=227)  1193    176        1 (1425/ 361-1571/ 438)   0.0 W08-W23 S48-S44
     16  626  263 (S30E55 L=156)  1342    155        1 ( 602/ 199- 649/ 350)  38.7 E64-E48 S35-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